Stable Fractional Vortices in the Cyclic States of Bose-Einstein Condensates 
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We propose methods to create fractional vortices in the cyclic state of an F — 2 spinor Bose- 
Einstein condensate by manipulating its internal spin structure using pulsed microwave and laser 
fields. The stability of such vortices is studied as a function of the rotation frequency of the confining 
harmonic trap both in pancake and cigar shaped condensates. We find a range of parameters for 
which the so-called jz -vortex state is energetically favorable. Such fractional vortices could be 
created in condensates of 87 Rb atoms using current experimental techniques facilitating probing of 
topological defects with non-Abelian statistics. 



The peculiar nature of quantum fluids stands out re- 
markably when observing their rotational characteristics. 
Rigid-body rotation is forbidden for a supcrfluid due to 
the constraint that its velocity field is bound to be ir- 
rotational. Instead they may acquire angular momen- 
tum by hosting a number of quantized vortex lines 
The superfluid state is characterized in terms of an order 
parameter which in general is a multicomponent com- 
plex function. In the single-component (scalar) case, the 
quantized vortex lines are phase singularities in the or- 
der parameter field, around which the complex phase of 
the order parameter changes by 2-kk, where k £ Z is the 
winding number of the vortex. 

In systems described by multicomponent order param- 
eters, the vortex lines may have more complicated struc- 
ture. A well-known example is the A-phase of 3 He, in 
which the Cooper paired fermionic atoms possess both 
orbital and spin angular momentum. The coreless half- 
integer vortex, for which the phase winding of the order 
parameter takes half-integer values, was discovered the- 
oretically in such system [2(. Existence of these kinds of 
topological defects is possible due to the discrete symme- 
try under combined gauge transformation and spin rota- 
tion of the related order parameter. Realization of Bose- 
Einstein condensates (BECs) in purely optical traps [3| 
has made it possible to create fractional vortex states 
also in ultracold atomic gases. In the absence of strong 
magnetic fields, the hyperfine spin degree of freedom of 
the atoms is unrestricted, leading to a multicomponent 
order parameter 

Condensates consisting of optically confined spin-2 
atoms are of special interest here. The first homotopy 
group 7Ti of the F = 2 manifold is non-Abelian [B| , en- 
abling non-commutable composition laws for two vortex 
lines [6( ■ Topological defects in the so-called cyclic state 
have been studied theoretically in [?J ■ The internal sym- 
metry of such order parameter can be mapped to the 
discrete tetrahedral group, enabling the existence of frac- 
tional 73 -vortices. In this paper, we describe two plau- 
sible schemes to create controllably a 73 -vortex into the 
cyclic state of an F = 2 spinor BEC. Moreover, we inves- 
tigate the energetic stability of such states in harmonic 



traps under rotation, paving the way towards realization 
of long-lived non-Abelian topological defects. 

An F = 2 spinor BEC is characterized by an order 
parameter which has five spin components, ipk, k — 
—2,-1,0,1,2, corresponding to the eigenstates of the 
five-dimensional representation of the spin operator S z . 
The different spin populations, Nk = j dr\ip k \ 2 : can be 
controlled accurately in an optical trap in the presence 
of an offset magnetic field @, Q. 

The ground state phase diagram of an F = 2 spinor 
BEC is divided into ferromagnetic, polar, and cyclic re- 
gions depending on the values of the spin-spin interaction 
strength [3 and the spin-singlet coupling constant 7 [lfj. 
This is in close analogy to the phase diagram predicted 
for d-wave paired Fermi systems Based on scatter- 
ing length measurements, the internal ground state of a 
87 Rb condensate belongs either to the polar or cyclic re- 
gion [13, Although several studies have predicted 
a negative value for 7, implying that the ground state 
would be polar, due to remaining experimental uncer- 
tainties it is not yet clear whether the ground state lies 
in the cyclic or polar region [1 31 ] - 

By restricting the atom population into the ip2 and 
if)— 1 components only, the value of 7 becomes insignifi- 
cant: The spin-singlet pairing interaction does not con- 
tribute to the energy since the spin-singlet pairing ampli- 
tude, 9 = Ylk(~ l) fe V'*V'-&3 vanishes. Hence the internal 
ground state of such condensate i s cy clic. In particular, 
let us consider states of the form [lj| 

where (r, ip, z) denote the cylindrical coordinates, and 
•d and 4> are functions of the azimuthal angle tp. 
Rzi^) — diag(exp[2i<^],exp[— i<j>}) describes a spin rota- 
tion through an angle <f> about the z-axis and e l1? is a 
gauge transformation. By choosing $ = <j) = ip/3, we 
obtain a state with unit (zero) phase winding in the 
■0+2 ("0-i) component. For such state, by traversing 
along a loop encircling the origin, the order parameter 
is continuously transformed onto itself accompanied by 
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a spin rotation of one-third of a full circle and a phase 
change of ^ . This configuration is therefore called a 73 - 
vortex state. Similarly, choosing -& — 2ip/3, <ft — — ^/3, 
yields a 2 / 3 -vortex state. 

Here we restrict to mean-field theory for which the 
zero-temperature energy of the condensate in the absence 
of external magnetic fields is given by (loj 



E = dr 



(2) 



with implied summation over repeated indices. In 
Eq. ©, h Q = -£-V 2 + F trap (r) - ft • L is the single- 
particle Hamiltonian where m is the mass of a 87 Rb 
atom, Vtrap the external optical potential, L the angu- 
lar momentum operator, and f2 the rotation frequency 
of the trap. The total density of the condensate is de- 
noted by n = iptipk, and the spin-spin interaction energy 
by (f) 2 = tp*'^%fijfli4>j'^u where is a matrix element 
of Ss/h. The parameters a, /3, and 7 are related to the 
s-wave scattering lengths in the total hyperfine state 
F according to a = ^(4ff2 + 354), (5 = — 7(52 — 94), and 

7 = l(9a ~ 54) - |(32 - where g F = ^-a F arc the 
bare coupling constants [lCj, [l5| . 

Minimizing the energy functional, Eq. @ , with a fixed 
total number of particles, leads to the zero-temperature 
Gross-Pitaevskii equation for the order parameter 

ha^k + an^ k + P(f)f 5 kl ^ + 7 6(-l) V- fe = Mo (3) 

where /i is the chemical potential of the system. We 
have solved Eq. ©, both in pancake and cigar-shaped 
harmonic traps, in order to find the ground states of the 
system. In the former case, we choose the axial trapping 
frequency w z large enough so that the condensate dynam- 
ics is frozen in the flattened direction, and hence the axial 
profile of the wave function can be accurately approxi- 
mated by a Gaussian function. In this limit, the interac- 
tion parameter a is replaced by the quasi-2D form a^D — 
2y/2ntiu! z a r [7(402 + 3<Z4)] . Setting the radial and axial 
trapping frequencies to (w r /27r, oj z /2tt) = (100, 1000) Hz, 
respectively, the harmonic oscillator length in the radial 
direction is a r = 1.1 /im, and we have chosen the value 
of a2D to correspond to 2 x 10 4 87 Rb atoms using the 
scattering length presented in [l(|. Moreover, we set 
(3 = a/100, (7 = —a/187) according to the values mea- 
sured in [131 ]. 

Figure Q] depicts the computed ground state profiles 
of the condensate for a range of rotation frequencies, 
fi, under the constraint that only the |+2) and | — 1) 
components have non-zero population. The top and 
bottom rows show the density of particles in the |+2) 
and | — 1) components, respectively. The field of view is 
20 a r x 20 a r , and the angular momentum of each compo- 
nent is given inside the figure. The horizontal bar below 
the figure shows the critical rotation frequencies for which 
the energies of neighboring ground state configurations 



are equal. The ground states are qualitatively similar for 
different values of f2 within each stability window, and 
the representative profiles shown refer to the rotation fre- 
quency at the middle of the interval. However, the total 
angular momentum and magnetization (f z ) of the con- 
densate varies within each stability window because the 
population ratio depends on the rotation frequency. For 
example, for the 73 -vortex state, illustrated in Fig. 1(b), 
the angular momentum ranges from 0.34 TVS to 0.40 NH 
and the magnetization from 0.031 N to 0.099 N. 

For the chosen parameters, the % -vortex state 
[ip + 2{r)e lLp , V-i ( r )] turns out to be energetically favor- 
able compared to the non-rotating ground state when 
fl > 0.0957 u r . As the rotation frequency is further 
increased above fl — 0.123 u r , a pair of 73- and %- 
vortices occupies the ground state. This state is analo- 
gous to the vortex molecule structure discussed in [l6| 
within the nonlinear sigma model. In the present con- 
text, the (Y3, 2 /s)-vortex pair can be viewed as a split 
gauge vortex, because the total spin rotation along a loop 
containing the pair vanishes and the gauge fields of the in- 
dividual vortices add up to an integer phase winding e lv . 
The 2 / 3 -vortex state carries roughly twice the angular 
momentum and kinetic energy of the 73 -vortex state. 
The former becomes energetically favorable compared to 
the non-rotating ground state at Q = 0.111 u) r , and at 
n = 0.127w r when compared to the y 3 -vortex state. 
However, this state is not the ground state of the sys- 
tem for any value of Q. For comparison, in the scalar 
\F = 2, = +2) condensate with the same parameter 
values, the singly quantized vortex state is the ground 
state within the interval fl € [0.164 u> r , 0.222 u> r ]. We 
have also confirmed that other possible vortex states in 
the given configuration, such as gauge or integer spin 
vortices, are energetically unstable. 

At even higher rotation frequencies additional 73 — 
and 2 / 3 -vortices enter the system, yielding a rich struc- 
ture in the ground state phase diagram, see Figs. QJd)- 
(h). With fast enough rotation, lattices consisting of 
and 2 / 3 -vortices appear. Vortex lattice transitions in 



such setup have been studied recently [17|. Moreover, 
vortex lattice structures in two-component condensates 
have previously been studied theoretically [HI, El, as well 
as observed experimentally [13] • In the present study, the 
population ratio is N +2 /N-i ~ 0.5, instead of being close 
to unity, which is inherited from the underlying mag- 
netic state (cyclic). For the parameter values of 87 Rb, 
we observe only interlaced square and rectangular vor- 
tex lattices. The shape of the unit cell depends on the 
rotation frequency as well as on the total particle num- 
ber N, tending more towards a square lattice in the limit 
Q — > uj ri and N — > 00. For example, with il = 0.90 u r 
and N = 2 x 10 4 , both unit cells are present: a square lat- 
tice close to the trap center where particle density is high, 
bounded by rectangular unit cells closer to the periphery 
of the cloud. The magnetization of the interlaced square 
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FIG. 1: (Color online) Ground state density profiles of 8r Rb atoms in the \F — 2, nip — +2) component (upper row) and 
\F = 2, m? = —1} component (lower row) in a harmonic trap rotating at frequency O under the constraint that only these 
components are populated. The density minima in the |+2) component correspond to fractional fa -vortices, and those in 
— 1) to 2 /3 -vortices. The l fa -vortex is the ground state within the interval Q £ [0. 0957 uj r , 0.123 uj r ] for the chosen parameter 
values. The color bar is scaled with respect to the maximum of the total density, max [n(r)]. 



lattice resembles notably that of the ground state of a 2D 
Ising model with antiferromagnetic interactions. More- 
over, these vortex lattices can be viewed as consisting of 
( V3 i % ) -vortex pairs, such as depicted in Fig. 0Jc). 

In the cigar-shaped limit, we choose the trapping fre- 
quencies (w r /27r, uj z /2tt) = (100, 10) Hz, and the coupling 
constants a, /?, (7) corresponding to 3 x 10 5 87 Rb atoms. 
The 73 -vortex state is the ground state within the inter- 
val SI € [0.171 Lu r , 0.237 o; r ]. The vortex line is bent in the 
stationary state of the system even in cylindrically sym- 
metric confinement, analogously to the scalar case [21|. 
Hence, the angular momentum of the 73 -vortex state 
is not determined only through the populations in the 
two components, N +2 and iV_i, but is strongly affected 
also by the bending of the vortex line. Within the sta- 
bility window, the atom number in the |+2) component 
ranges from 0.33 N to 0.41 N, whereas the total angular 
momentum ranges from 0.21 TVS, to 0.39 NH. For com- 
parison, the singly quantized vortex state in the scalar 
\F = 2.niF = +2) case is energetically favorable within 
the interval SI € [0.294 w r , 0.394 u r ]. 

As the rotation frequency is increased further, the 
ground state contains a ( 1 /z , 73 ) -vortex pair which is il- 
lustrated in the isosurface plots of Fig. [2]for SI = 0.25 u> r . 
At each point on the surface shown in Fig.^a), the par- 
ticle density in the |+2) component is equal to 10% of 
the total particle density, and similarly for the | — 1) com- 
ponent in Fig. [2f b) . The surfaces are cut open to clarify 
the vortex lines inside the clouds. 

The ground state energy as a function of rotation fre- 
quency SI is plotted in Fig. f^c) (the solid curve). The 
interval within which the 1 fa -vortex is stable is bounded 
by the dashed vertical lines. The slope of the energy 
curve steepens as the total angular momentum of the 
state increases. The angular momentum of each com- 
ponent is shown with circles for |+2) and triangles for 
| — 1). Angular momentum of the | — 1) component devi- 



ates slightly from zero even in the 73 -vortex state since 
bending of the vortex line breaks the cylindrical symme- 
try in both components. For higher rotation frequencies, 
a square lattice of vortices forms also in the cigar-shaped 
case. 

Finally, we will describe two alternative methods of 
creating a 73 -vortex state in a trapped condensate. The 
first scheme is a direct generalization of the method uti- 
lized by Matthews et al. to create the first vortex in 
a two-component BEC 2^, 23|. The condensate is first 
prepared in the \F = \,mp = 0) state in an optical 
trap, e.g., by applying a magnetic field gradient during 
the evaporative cooling stage [23|. Subsequently, two- 
photon microwave and auxiliary laser fields are used, as 
in [23}, to transfer roughly one-third of the population 
directly into a finite orbital angular momentum state in 
\F = 2,iriF = +2). Thereafter, the desired configura- 
tion in the F = 2 manifold is achieved by switching off 
the laser beam and transferring the remaining popula- 
tion from \F = I, nip = 0) into \F = 2, nip = —1) by 
a resonant microwave field. The scheme described above 
yields a l fa -vortex state [ip +2 (r, z)e lip , tp_i (r, z)] with the 
desired population ratio of N + 2/N-i ~ 0.5. Following 
the same prescription but initially preparing the atoms 
in \F = l,mp = —1) results in [ip + i(r, z)e zip , "0-2(7, z)\ 
which is a 73 -vortex state. However, in this case some 
of the initially prepared atoms would have to be dis- 
carded in order to obtain the optimal population ratio 
with maximal transfer of 50% with the two-photon mi- 
crowave field. The auxiliary far-detuned rotating laser 
beam could in principle be utilized in rotating the con- 
densate [111, providing stability for the fractional vortex 
state. 

Another possibility to create the 73 -vortex state em- 
ploys a Laguerre-Gaussian (LG) laser mode which has 
previously been utilized for creating multiply quantized 
vortices in a scalar BEC as well as vortex- antivortex su- 
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FIG. 2: (Color online) (a) Isosurface plot of the density of 
particles in the |+2) component |i/>+2(r)| 2 , and (b) in the | — 1) 
component |i/>_i(r)| 2 . The field of view is illustrated approxi- 
mately by the double arrows. In both isosurfaces, the density 
is equal to 10% of the maximum total density of particles, 
max [n(r)]. For the chosen rotation frequency of Q — 0.25 u) r , 
the ground state contains a (V3 , % ) -vortex pair. The vor- 
tex lines, for clarity drawn with different color as the rest of 
the isosurface, are bent even in an axisymmetric trap, (c) 
Angular momentum in the |+2) component (circles) and in 
the | — 1) component (triangles) as a function of the rotation 
frequency Q. The solid curve shows the ground state energy 
of the condensate. The 73 -vortex is the ground state config- 
uration between the dashed vertical lines. 



perposition states in spinor condensates [26J, l27(. The 
73 -vortex configuration could be achieved by preparing 
the condensate initially in \F = l,mp = 0), and trans- 
ferring roughly one-third of the atoms into \F = 2, mp — 
+2) using counterpropagating a + polarized LG and a~ 
polarized Gaussian laser beams [271 • The remaining pop- 
ulation in \F — l,mp — 0) could be transferred into 
\F = 2,mp = —1) using resonant microwave pulses yield- 
ing the desired vortex state and population ratio. 

In experiments, the lifetime of a condensate in the un- 
stable F = 2 manifold has been limited to a few hundred 
milliseconds due to inelastic collisions Therefore, it 
seems more promising to create the 73 -vortex dynami- 
cally instead of relying on the cloud to thermalize into the 
lowest energy state. The magnetization relaxation time 
scale is beyond current lifetimes of the condensates, and 
hence the population ratio of the |+2) and | — 1) compo- 
nents may be expected to be conserved even in the pres- 
ence of weak magnetic stray fields. Moreover, very slow 
spin-mixing dynamics has been reported in 87 Rb conden- 
sates initially prepared in the cyclic configuration [8|, in 
favor of stable fractional vortex states. 

In summary, we have proposed two realistic meth- 



ods to create a fractional 73 -vortex state in an F = 2 
spinor Bose-Einstein condensate in a laboratory experi- 
ment. Both schemes are based on transferring popula- 
tion from the stable F = 1 manifold into two spin states, 
\mp = +2) and \mp = —1), in the F = 2 manifold by 
using a combination of laser beams and pulsed magnetic 
fields. We have performed energetic stability analysis of 
such fractional vortex states, and found that, e.g., the 
73 -vortex is the energetically favored state within an in- 
terval of rotation frequencies of the confining potential. 
For increasing rotation frequencies, the ground state con- 
tains more and more of 73 - and 73 -vortices, which even- 
tually form an interlaced square lattice. The energetic 
stability of the fractional vortices under rotation raises 
hope for experimental realization of topological defects 
with non-Abelian statistics. Effects related to the non- 
Abelian nature of the fractional vortices could be probed 
for instance by vortex- vortex collision experiments @ . 
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